Structural properties and quasiparticle energies of cubic SrO, MgO and SrTiOa 



Giancarlo Cappellini 

INFM Sezione di Cagliari and Dipartimento di Fisica, Universitd di Cagliari, Cittadella Universitaria, Strada Prov.le 

Monserrato-Sestu km.0.700, 1-09042 Monserrato (Ca), Italy 

Sophie Bouette-Russo, Bernard Amadon, Claudine Noguera and Fabio Finocchi 
Laboratoire de Physique des Solides, UMR CNRS 8502, Batimemt 510, Universite Paris-Sud, 91405 Orsay, France 



The structural properties and the band structures of the charge-transfer insulating oxides SrO, 
MgO and SrTiOa are computed both within density functional theory in the local density approxi- 
mation (LDA) and in the Hedin's GW scheme for self-energy corrections, by using a model dielectric 
function, which approximately includes local field and dynamical effects. The deep valence states 
are shifted to higher binding energies by the GW method, in very good agreement with photoemis- 
sion spectra. Since in all these oxides the direct gaps at high-symmetry points of the Brillouin zone 
may be strongly sensitive to the actual value of the lattice parameter a, already at the LDA level, 
self-energy corrections are computed both at the theoretical and the experimental a. For MgO and 
SrO, the values of the transition energies between the valence and the conduction bands is improved 
by GW corrections, while for SrTiOs they are overestimated. The results are discussed in relation 
to the importance of local field effects and to the nature of the electronic states in these insulating 
oxides. 
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I. INTRODUCTION 



Aside from being the major constituent of the outer crust of the earth, oxides are of considerable interest for their 
wide applications in various technological fields such as electrochemistry, catalysis or microelectronics!!!. They present 
very diverse electronic ground states which range from superconducting, metallic, semi-conducting (intrinsic, charge 
density- or spin density-like) to insulating. ^From a theoretical point of view, a proper description of their electronic 
properties is still an active research subject. On the one hand, several transition metal oxides require an inclusion 
of many-body effects beyond the available effective one-electron approaches, such as the density functional theory 
(DFT) or the ab initio Hartree-Fock (HF) method. On the other hand, even for insulating charge transfer oxides, 
such as MgO and SrO, whose electronic ground state is well accounted for in a one-electron picture, the calculation 
of quasiparticle (QP) energies within the DFT or the HF method gives rather poor results. Indeed, one has to face 
the so-called energy band gap problemEl. Indeed, when using the Kohn-Sham eigenvalues of the highest occupied and 
the lowest empty states, the resulting gap usually underestimates the experimental data. On the other hand, the 
corresponding difference between the HF eigenvalues overestimates the gap severely. Alternative ways to calculate 
the band gap in infinite systems have been found, in the framework of HF or DFT. One of them is based on total 
energy differences between ground and charge transfer statesa. [The other one relies on the value of the gap found in 
the one electron spectrum of the ionized or hole-doped systemcl. However, such methods are unable to provide the 
full QP spectrum. 

In the last decade, due to the implementation of angular resolved photocinission experiments able to cope with 
charging effects on insulating oxides, detailed QP spectra of MgOQ, SrTi03@ or Ti02&B among other oxides, have 
become available. The experimental achievements demand systematic improvements in the theoretical description of 
the electronic excitations of simple oxides. 

In semiconducting and ionic crystalline compounds, corrections to the DFT eigenvalues, as obtained from the Kohn- 
Sham equations, may be estimated by first-order many-body perturbation theoryBtJ, with respect to the difference 
between the exchange-correlation (XC) self-energy and the corresponding XC potential of the DFT, for which various 
approximations arej-available. Within the same approach, the self-energy operator is treated in the GW approximation 
proposed, by, HedinO. Such calculations have been carried out with success in several semiconducting and insulating 
crystalai3il3, permitting the prediction of the quasiparticle energies and the calculation of optical properties of real 
solids. 

One of the crucial points in such calculations is the evaluation of the self-energy operator, which is numerically very 
expensive, because of the required knowledge of the full inverse dielectric matrix of the system under study. This 
is a bottLeiieck for the application of the method to systems with a large number of atoms in the unit cell, such as 
surfaccso'EJ and defective or low-symmetry solids. In order to reduce the remarkable numerical effort and be able 
to predict electronic excited states in a wider range of real systems, one may use efficient GW methods, in which a 
model dielectric function is used to mimic the screening properties of the system under study. I— . 

Pioneering work in this direction has been performed by Bechstedt and Del Sole in a tight-binding schemetS 
and by Gygi and Baldereschilj within the DFT in the Local Density Approximation (LDA) for the exchange and 
correlation energy. The key of their work relies (i) on the use of a model dielectric function, which mimics the 
main screening properties of semiconductors, and (ii) on an approximated self-energy operator in which both local 
fields and dynamical screening - giving rise to opposite contributions to band gaps- are neglected.E3HL3 For a large 
family of semiconductors, the band gaps obtained through these simplified methods show good accordance both with 
experiments and full GW results, in which the screening properties of the system are calculated from first principles. 
More recently, an efficient GW method that approximately-includes dynamical-screening and local field effects, without 
increasing the computational effort, has been formulatedE3 and proved to reproduce quasipapticle energies and band 
gapsJbr either bulk zmpblende solids (ranging from Si, GaAs, AlAs to the more ionic ZnSe)ci or non cubic systems 
(BNo, SiC polytypesoo) satisfactorily. The method is based on a linear expansion of the dynamical contribution to 
the exchange-correlation self-energy, on the use of a model dielectric function, and on a local ansatz for the treatment 
of the intrinsically non local screened interaction. Depending on the system, the speedup of the calculation may result 
of two order of magnitude with respect to corresponding full GW ones. Moreover, this kind of GW calculations is 
feasible on a standard workstation. 

The aim of the present paper is to study the ground state properties of cubic magnesium oxide, strontium oxide 
and strontium titanate, in the DFT-LDA scheme, as well as their quasiparticle band structures within the efficient 
GW scheme. On one hand, no previous GW calculation of SrO or SrTi03 exists, to our knowledge. On the other 
hand, these charge transfer oxides are characterized by different fundamental gaps, ranging from ~ 3 eV (SrTiOa) 
to ~ 7.7 eV (MgO), as given by experimental data. Their bandstructures may substantially differ from those of 
the semiconductors cited above, due to the presence of strongly localized electronic states. Thus, such a study may 
also represent a test for the approximate GW scheme, originally proposed for small-gap semiconductors, to show its 
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capabilities relatively to a different class of physical systems. 

The paper is organized as follows. The theoretical framework is depicted in Section II and it is followed, in Section 
III, by a description of the computational ingredients used for the description of structural and electronic properties 
of MgO, SrO and SrTiC>3. The results obtained for these systems are presented separately in Sections IV, V and VI, 
respectively. Finally, in Section VII we present our conclusions and the perspectives opened by our work. 



II. GW CORRECTIONS 

The energy J2 £ of a quasiparticle in a band n at a given k point in the Brillouin zone can be obtained through the 
equationM^^! 



faW + ^W,r';E nt )^{r') = E n ^ k {f) (1) 



which is written in atomic units (e=Ti =m e =l). V ey ±(r) is the external potential originated by the ions, Vn(r) is the 
Hartree potential and X is the mass operator, usually called the self-energy operator, which is in general non local, 
non-Hcrmitian and energy dependent. In principle it contains the exchange rand correlation effects and is state (nk) 
dependent. One can compare the above equation to that of Kohn and Shamed: 



-^- + V cxt (r) + V li (r) + V xc (r) < ~> - 



*nSW = E nl^l (2) 



where V xc is the exchange-correlation potential of the DFT (for which several approximation are available) , and the 
superscript will denote the quantities computed at the DFT level from now on. E^ } is usually considered as a first 
approximation to the quasiparticle energy. The failure of this scheme for evaluati ng tran sition energies E ,? — E t 
between different quasiparticle states in solids has been evidenced by many authors9H30. 

One can correct the DFT eigenvalues E g, in a first-order perturbation theory with respect to £ n g — V xc , and 
estimate the QP energies as: 

E nk * E< 2+ < ^nk\^' E nk) ~ MOI^g > ( 3 ) 

The evaluation of S(r, r'\ E A is in general a very difficult task. A possible approximation is the GW one,i in which 
a perturbation expansion for the self-energy can be constructed and stopped at the first order: 

/+OO J I 
e +i8u G{r,?;u + u')W{r,r';u') (4) 
-oo 2?r 

where G is the one- particle Green function, W is the screened Coulomb interaction and 5 = + . Following Bech- 
stedt and coworkers ,&El the real part of the self-energy operator (which is relevant for the corrections to the QP 
bandstructure) may be separated in a dynamic and a static contribution, and written as: 

S E(f, r» = E scx (f, r') + S coll (f, P) + £ dy,1 (r, r>) (5) 

S scx and S coh are the static screened-exchange (SEX) and the coulomb-hole (COH) terms, respectively, and E dyn the 
energy dependent contribution. Although the previous equation is purely formal, it is the starting point of simplified 
GW schemes, since one can work out approximated expressions for each contribution separately. The most expensive 
step towards the evaluation of GW corrections to £l' i is the calculation of the energy-dependent screened interaction 
W(r, r';u>). Although many i mpr ovements have been proposed for an efficient calculation of the dielectric response 
function from first principles ,EZrc3 still this stage is both computer time and memory very demanding when a large 
basis set is used for the description of the electronic structure. In particular, this is the case of oxides in a plane- wave 
pseudopotential approach. On the other hand, the improvement of DFT-LDA band structures obtained by using 
schemes originally formulated on phenomenological basis, such as the scissor operator or the Slater's a potential,EI 
showed that the essential physics can be described through the use of effective or state-dependent potentials. From 
a more fundamental point of view, one can work out a model for the screening properties of the system, and look to 
what extent the computed QP spectra account for the observed experimental values. Possibly, the mo del .should also 
be simple enough to be generalized to a large number of different systems. Following previous studies,EjH we use a 
model dielectric function to estimate both the static and the dynamical contributions to the self-energy (Eq.^), since 
it is the key ingredient entering the efficient GW approximations. 
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A. Static Screening and Local Field Effects 



For trie-diagonal part of the screening function, we use an analytical model originally proposed by Bechstedt and 
Del SolerEa 



e(q,p) = 1 + 



3<Z 4 



1 



Itf 



(6) 



where the Fermi vector kp and the Thomas-Fermi wave- vector qxF depend on the average electron density p. This 
expression mimics the free electron gas behavior at high q. At q = 0, Eq.|| gives the value of the optical dielectric 
constant Eoo of the oxide. Moreover we use a local ansatz for the static Coulomb screened interaction, analogously to 
that of Hybertsen and LouieE3: 



W(f,P) = \[W h {f-f;p{r)) + W h (f-r';p(f'))} 



(7) 



where W h is the screened Coulomb interaction of a virtual homogeneous system characterized by a finite optical 
dielectric constant. Ei£H_Duasiparticlc energies determined from the above ansatz result in good accordance with full 
GW calculated ones.otZl It permits, together with Eq.(6), to obtain a simple analytical-expression for the static part 
of the Coulomb hole self energy S coh . The explicit formula has been given elsewhere.Ell For the calculation of E sex 
in Eq.[|, the Fourier transform of W(r, f') is needed, namely W(k + G,k + G'). In order to reduce sharply the time 
needed to compute this term, we proceed in two steps. Firstly, we retain only the diagonal terms in the screened 
coulomb interaction. Secondly, we account for local field effects (i.e. non diagonal terms with G ^ G' in W) by using 
suitable state-dependent densities to compute kp and qTF- Given a state nk, we define the effective electron density 
which determines the typical screening lengths q^p and kp 1 , relevant for the estimation of E g (Eq.^) as the 
expectation value of the ground state electronic density p(r) on the actual state: 



d 3 rp(r)\ip nj :(r)f 



(8) 



Since the relevant screening lengths kp 1 and q^p are increasing functions of p n ^, the latter is intended to give an 
approximate description of the local fields in E sox through an enhanced screening for the electronic states which 
contribute effectively to the total electron density. In the oxides considered here, the effective densities p % range 
from few tenths of the mean valence electron density p = Z va i/V for conduction states up to 4-5 p for deep valence 
states. As a consequence, the lengths kp 1 and q^p are shorter for the screening of valence band states than for 
conduction band states. This difference is enhanced with respect to more covalent semiconductors, such as Si and 
GaAscJ, in which the valence electron distribution is more delocalized and less dissimilar from that of the conduction 
states than in oxides such as MgO. 

This approximate scheme for treating local field effects was analyzed in detail and its reliability was. shown through 
a comparison with full GW calculations, in the case of both small- and wide-gap semiconductorscJ'El. In Ref. |2l], 
it was also shown that the use of a non-diagonal model dielectric function does not clearly improve band gaps while 
bringing additional computational costs. 



B. Dynamical contributions 

The diagonal matrix element of the dynamical part of the self energy (the third operator on the right side in Eqj^) 
on the state ip p 

S S n H = ^|S dyn (f,P;c)|^) (9) 
can be linearly expanded in energy around the DFT value E } , with linear coefficient: 



This approximation for the energy dependence of the self-energy xuperator turns out to be valid in a wide energy range, 
as confirmed by the results of full GW calculations on bulk Si.E3 Solving the Dyson equation within the first-order 
perturbation theory with respect to S — V KC , and using Eqs.|| and [|, the QP shifts for the ip n k state hence read: 
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S c<A + E sex + E dyn,£(0)x _ \y \ 

u * \t riK i ' n/c ' 



771 — — — n ^ nfc xc nk' fll) 

nk nk - 1 + r 

^nk 



where ^ co i l and £ sc -? are the expectation values of the static COH and SEX contributions to the self energy, respectively. 
^From the GW expression of £, replacing G by G^ **, computed at the DFT leveled : 

tl k n, k nk tl k 



n'k' 



where 



£»1 = / ^r^(r)e^-V n , fc -,(r) (14) 



The model e 1 is diagonal in g space and the local field effects are accounted for through the use of the state averaged 
density p n ^ (Eq.^|). The (^dependent dielectric function is obtained from the static one by means of the plasmon- 
pole approximation (PPA)ll3. The PPA has been widely used to compute semiconductor and insulator quasiparticle 
bandstructures, even in the energy range of semicore statesSH In order to account for the differences E^°J — E^L 

OJ ° nk n'k' 

in Eqsjl2| and |i~3| , replacing them with the state independent form E ^ — E , L w —q 2 /2 as in Ref. El and using the 
sum rule J2 n 'k' \ B nn'(^)\ 2 = 1j one obtains: 

»* v n fe ; y (27r )3 y wvr w + £ ^ 

P ^-J (2^)3 g 2 „ (w + £ )2 w 

Therefore, once the spectrum and eigenfunctions obtained in the framework of the DFT are known, the calculation 
of E d g" (E^ )) and @ nk - according to Eqs.[l|| and needs no additional external parameters apart from the model 
dielectric function, and results a simple task that can be carried out on a standard workstation. 

We have made extensive checks to test the reliability of the approximations used to obtain Eqs.[l5] and |l6|, in the 
case of Silicon. Some transitions between high-symmetry points of some valence and conduction band states in Si 
are compared in Tables I to those obtained through full GW. The energy of the transitions between high-symmetry 
points, as well as the fundamental gap in Si, are in very close agreement with those obtained from full GW calculations 
and with experiments for both methods. As far as the QP eigenvalues in Si are concerned, the top of the valence 
band r25',„ is shifted downwards with respect to the DFT values of 0.28 eV. This can be compared with the upward 
shift of 0.07 eV, as found by Godby et alH. 



Si 


full GW (a) 


full GW (b) 


this work 


Exp. 


— MSB', v 


12.04 




12.09 


12.5±0.6 


^25',v — Tis^c 


3.35 


3.30 


3.33 


3.4 


rw, v — IV, c 


4.08 


4.27 


4.25 


4.2 


Ij2' ,v — ^25' ,v 


9.79 




9.98 


9.3±0.4 


L\ >v — ^25', v 


7.18 




7.27 


6.7±0.2 


L>3' ,v ~ ^25' ,v 


1.27 




1.25 


1.2±0.2, 1.5 


T*25',v — Ll,c 


2.27 


2.30 


2.22 


2.4±0.2, 2.1 


^25',v — ^3,c 


4.24 


4.11 


4.15 


4.15±0.1 


L3 1 ,v — il.c 


3.54 




3.47 


3.45 


L3' ,v — L3,c 


5.51 




5.41 


5.50 


Xi lV — T25' :V 


2.99 




2.98 


3.3±0.2, 2.9 
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1.44 




1.54 


1.3 


v — c minimal gap 


1.29 


1.24 


1.21 


1.17 



TABLE I. Computed vertical transitions and minimal valence-conduction band gap for Silicon (in eV). (a): Ref. [27j (b): 
Ref. [111. The experimental values are quoted in Ref. E71 
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III. COMPUTATIONAL DETAILS 



The Density Functional calculations are carried out within the LDA for exchange and correlation!^. The Kohn- 
Sham orbitals are expanded in a plane-wave basis set. Special care has bees, used in constructing the pseudopotentials 
for cations (Sr, Ti, Mg), in order to avoid the occurrence of ghost statesc3 and to assure an optimal transferability 
over a wide energy range. We find that the inclusion of semicore states, such as 4s and Ap states for Sr, and 3s and 3p 
states for Ti, greatly improves the transferability of the pseudopotentials, and is at the same time necessary to take 
into account the hybridization of cation semicore states with O 2s states. .— . 

Angular components up to I = 2 are included. The scheme of Martins and Troullier is usedEJ to generate separable 
soft norm-conserving pseudopotentials, with core radii (bohr): 2.00 (Sr, 4s), 1.50 (Sr, Ap), 1.90 (Sr, 4d); 1.30 (Ti, 3s), 
1.40 (Ti, 3p), 1.80 (Ti, 3d); 1.38 (O, 2s), 1.60 (O, 2p), 1.38 (O, 3d); 2.00 (Mg, 3s, 3p, 3d). The cutoff energy needed 
to obtain a convergence better than 0.1 eV of both total energy and Kohn-Sham eigenvalues is found to be equal to 
50 Ry for MgO, 60 Ry for SrO aniL80 Ry for SrTi0 3 . 

The use of 10 special k pointsc.3 for charge integration in the irreducible wedge of the Brillouin Zone (IBZ) is 
sufficient to achieve a good accuracy for the computed total energy; for instance, the total energy changes by less 
than 10 meV and the fundamental band gap by less than 0.1 eV when using 10 instead of 6 special points to sample 
the IBZ. 

The main bottleneck in the ab-initio GW calculations is the calculation and the inversion of the full dielectric matrix 
from the eigenfunctions and eigenvalues of the system. Typically, this task takes approximatively 75% of CPU time 
needed for the calculation of the self-energy correction for a single statea. As a result of the use of a model screening, 
the computational cost and memory requirement are strongly reduced, and are comparable to those of calculations 
carried out within the DFT. Efficient GW calculations can thus be carried out also on a common workstation.!^ 

The evaluation of the GW corrections to the LDA bandstructures, and especially of the SEX contribution S scx to 
the self-energy, needs some care, due to the presence of an integrable divergence. A reduction of the numerical cffotjb 
by using a limited number of k points in the IBZ is possible through the method proposed by Gygi and BaldereschiE3. 
In our calculations, the regularization is performed, for all the fc-points for which £ SGX is calculated, by transforming 
the integral over the BZ of the diverging contribution into an integral of a function F(k) (periodic over the BZ), 
that can be computed analytically, plus a discrete sum over special fc-points of a smooth function. For the sake of 
conciseness, in the case of the bare exchange term (the generalization to the case of the diagonal screened exchange 
in cubic systems is straigthforward), the transformation reads: 

Jbz (2t0^ \k'~k~g\ 2 ~ ^ J bz W ^ } ^"'^ + 4^' 

1 al k' 

(17) 

9 BZ (q) = 1 if q € BZ and zero otherwise, B*£,(g) has previously been defined (EqJT4[), and uig is the weight 
associated to the special point k in the IBZ. By following this method, 10 special k points in the IBZ are sufficient to 
get converged E sex within few tens of meV. 

Because of the presence of the static screened exchange term, each GW correction to a given E^°), scales as NrNNp™, 

nk' K 

where Nr, N and iVpw are the number of k points used in the summation of the charge, of occupied electronic states, 
and of the plane waves used in the expansion of the Kohn-Sham orbitals, respectively. 



l ?r'r^i - o BZ (k'-k-g) \BZm 2 F(k'-k) 

k' — k — q\ z 



IV. MAGNESIUM OXIDE 

The ground state properties obtained for the 0\ phase (rocksalt) of magnesium oxide are summarized in Tab.[n|. 
They are computed by fitting the curve of total energy versus lattice parameter to the equation. of | st ate proposed by 
Murnaghan. One can note that our results compare well to other all-electron LDA calculationsllSEZl and show slight 
discrepancies with respect to experimental data: the computed lattice parameter is underestimated by about 2% and 
the cohesive energy is about 15 % too large. We attribute these errors mainly to the use of the LDA. 

The bandstructures, both in the LDA and including quasiparticle corrections, are shown in Fig.|l| along high- 
symmetry lines in the BZ. In the following, we refer to the energy levels relative to the top of the valence band in 
the LDA calculation, which is arbitrarily set to zero. The bands are calculated at the theoretical equilibrium lattice 
parameter ath=4.125 A. In Table III we detail the energy differences between electronic states at high-symmetry 



points in the BZ, computed either at the theoretical equilibrium lattice parameter ath=4.125 A or at the experimental 
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one a OXD =4.211 A. One can easily see that some of the direct transitions, and in particular those at T, strongly depend 
upon the value of the lattice parameter. This is true at the LDA level, while the GW corrections are less sensitive to 
the actual value of a. In ionic rocksalt compounds, the value of the fundamental gap at the T point is indeed driven 
by the strong Madelung potential, which varies as the-ipverse of the lattice parameter. 

In agreement with a previous full GW calculationlia, the quasiparticle corrections consist of a rigid shift of the 
conduction band (CB) with respect to the valence band (VB) — that is, application of a scissor operator, which is 
modulated by smooth k dependent terms of the order of 10 % of the rigid shift itself. The total valence bandwidth 
ri5,u — Ti^ changes from 17.3 eV in LDA to 20.1 eV with quasiparticle corrections. Similarly, the V and X\ points 
move ~ 2 eV downwards from the the valence band maximum (VBM), as one can see from Fig.Q. This improves the 
accordance with the experimental XPS spectra which show a main peak at 18 eV and a shoulder at 21 eV below the 
VBMEl The fundamental band gap passes from 5.21 eV to 8.88 eV. .— ■ 

In all these calculations, we assume the experimental value = 2.95cJ for the static dielectric constant. Recently, 
through an ab initio calculation of the dielectric function within the LDA, Shirley found = 3.03, at the experimental 
lattice parameter^]. In order to test the sensitivity of our theoretical scheme to the choice of eoo, which enters as a 
parameter in our model screening function defined in Eq|| we performed an additional GW calculation by adopting 
£oo = 3.50. Consistently with a stronger screening, the direct transition energies at high-symmetry points are lowered 
by ~ 0.5 eV at most. 

It is interesting to note that the main corrections to the Kohn-Sham eigenvalues of VB states come from the 
static screened exchange term E sox (Eq.||), while the coulomb hole contribution £ coh dominates the quasiparticle 
corrections to the eigenvalues of the CB. For instance, Sp™ ~ — 16eV and E r ° h ~ —10 eV at the bottom of the VB, 

while E?f x ~ — 5 eV and £ r oh ~ — 8 eV at the bottom of the CB. This is consistent with the more localized nature 
1 1 1 1 

of the VB states with respect to the CB states. The same remark applies to SrO and SrTiOa, which are discussed in 
the following sections. 

As far as the interpretation of the main optical transitions is concerned, we indicate-in Table III the experimental 
values with their tentative assignement, according to Schomberger and AryasetiawanEj. We agree with their inter- 
pretation, since our GW values, computed at the experimental lattice constant, differ by ~ leV at most from theirs. 
Our calculations slightly overestimate the experimental transition energies. 



MgO 


a(A) 


B(Mhar) 


E coh (eV) 


Experiment 


4.211 


1.55^ 


10.33 


Ref. 
Ref. 
Ref. 


42 

if 

37 




4.191 

4.16 

4.09 


1.46 
1.71 


10.67 


Present work 


4.125 


1.56 


11.80 



TABLK-JI. Structural prpmefjties of Cubic MgO. A comparison is made with other LDA calculations, in a pseudopotential 
approaches or using LMTOtlO 



MgO 




£ g (r) 


E S (X) 


E g (L) 


Minimal Gap 


Experiment^ 


7.7 


13.3 


10.8 


7.7 (r - r) 


HF+p«lJd 

ldab n 

LDArfjGWO 

ldaH n 

LDA+GWC2 


8.21 

5.2 

7.7 

4.73 

7.81 


15.79 

10.5 

13 


11.48 

8.9 

11.4 


8.21 (r - r) 
5.2 (r - r) 
7.7 (r - r) 
4.73 (r - r) 
7.8i (r - r) 


this work 

LDA 1 " 

LDA* 

LDA+GW* 
LDA+GW* 


5.21 
4.61 
8.88 
8.20 


10.40 
10.26 
14.43 
14.22 


8.93 
8.46 
12.99 
12.46 


5.21 (r - r) 
4.6i (r - r) 
8.88 (r - r) 
8.20 (r - r) 



TABLE III. Direct gaps of Cubic MgO (in eV). Our results are given both in the LDA and with quasiparticle corrections 
(LDA+GW), either at the theoretical equilibrium lattice constant (') or at the experimental one ('). The theoretical calculations 
of Refs. |44[|l5|,|l0| are all carried out at the experimental lattice constant. 
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V. STRONTIUM OXIDE 



The ground state properties computed for the phase of SrO are summarized in Tab.IV. The bandstructures, 
computed in LDA and with QP corrections, are shown in Fig.^| along high symmetry directions. 

A deep band, originating mainly from Sr 4s states, is not shown in the figure. It is found around -30.5 eV in the 
LDA calculation, with a dispersion less than leV. The QP corrections shift this band to lower energies, about 33 eV 
below the valence band maximum (VBM). 

The bands of lowest energy shown in Fig.|| originate from Sr 4p and O 2s states. In the LDA, the bandwidth 
is about 4.3 eV, while the GW corrections push that band towards higher binding energies, and split it in two less 
dispersive structures, of which the lower originates from O 2s states mainly, and the upper from Sr Ap states. The 
splitting between O 2s and Sr Ap levels is essentially due to the larger E sex contribution for the former states, as a 
consequence of their stronger localization. i_ . 

As far as the deep VB levels are concerned, XPS experimental found two broad peaks at ~ —35 eV (A) and ~ 
-17.5 eV (B) below the top of the valence band. The authors assigned the A structure to the band arising from Sr 4s 
states, and the B structure to those coming from O 2s and Sr 4p states. The resolution, however, was too poor for 
them to resolve the two latter contributions to peak B. The positions of the deep VB levels are badly accounted for 
in the LDA, which would yield the A peak at -30|.5 eV and a broader B structure around -13 eV. On the other hand, 
Hartree-Fock calculations including correlation^ give a separation of about 6.5 eV between the O 2s and the Sr 4p 
states, which is likely overestimated. A better agreement with experimental data is obtained through the inclusion of 
GW corrections, which shift the Sr 4s band to ~ -33 eV and give two contributions to the electronic density of states 
at ~ -18 eV and ~ -14 ey-,with respect to the VBM, about 0.5 eV wide. The latters seem consistent with the broad 
B structure seen in XPSE3. A later XPS investigations found a splitting Aab between the A and the B structure 
equal to 17.9 eV, which was interpreted as the difference in the binding energy of Sr 4s and Ap levels in SrO. Since it 
seems difficult to disentangle the Sr Ap and the O 2s contribution to the B structure, we compare the minimum and 
the maximum splitting Aab as obtained in our calculations. In the LDA, Aab ranges from 14.5 eV to 17.5 eV, while 
that obtained by including the GW corrections ranges from 16 eV to 19 eV and better agrees with the XPS data. 
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The upper part of the valence band results very similar in the LDA and GW calculations. Its width is 2.2 eV 
in LDA and 2.1 eV after GW corrections. The dispersions of both QP and LDA eigenvalues along high-symmetry 
directions are very similar, too. 

The lower part of the conduction band drawn in the figures is shifted upward by GW corrections, opening the LDA 
gap by about 3 eV, with additional and small corrections depending on the particular band and k point as one can 
see in Table @. In particular, the direct gap at P (X) passes from E S (T) = 4.29 eV (E g (X) = 3.03 eV) in LDA to 
-Eg(P) = 7.54 eV (E S (X) = 6.39 eV) after GW corrections. It is important to note that the direct gap P15 - Li is 
strongly sensitive to the precise value of the lattice parameter ao: while the above value E S (T) — 4.29 eV is calculated 
at the LDA equilibrium lattice parameter ath, at a C xp = 5.16A we find a smaller E g (T) — 3.92 eV, i.e. a reduction of 
about 0.4 eV. The behaviour at the X point is different, since the E g {X) calculated at a CX p is equal to 3.11 eV, about 
0.1 eV bigger than that computed at ath- Given the different nature of the electronic states at the various symmetry 
points in the BZ, a non trivial dependence of the eigenvalues upon the structural parameters is generally expectecO. 

On the other hand, the QP energy differences obtained by our method are less dependent on the precise value of the 
optical dielectric constant used in the model screening (Eq.^). For instance, E g (T) varies only by -3% when passing 
from = 3.35 to £oo = 3.7, which represents a 11% increase of the optical dielectric constant. 

We are not aware of any angle-resolved photoemission or inverse photoemission experiments on SrO, able to provide 
direct experimental information on the band dispersion around the Fermi level. Thus, the comparison of our resulte 
to experimental data is mainly based on the information issued from reflectivity measurements. Rao and KearneycZI 
conclude for direct gaps at T (X) E g (T) = 5.9 cV (E g (X) — 6.28 eV), and exciton binding energies at high-symmetry 
points of the RZ around 0.2 - 0.3 eV. Conversely, more recent optical measurements on well characterized single 
crystals of SrOc3 found E g (X) = 5.79 eV, lower than E g (T) = 6.08 eV. On the basis of the onset of the optical 
absorption spectra the authors also deduce that SrO has an indirect minimum gap, while BaO shows a direct gap. 
at X. On the .theoretical side, the nature of the fundamental gap of SrO is still debated: while both LMTO/LDAHI 
and APW/XqEJ calculations give an indirect gap Ti 5 — X 3 clearly underestimated (3.8 and 3.9 eV respectively), 
Pandey and coworkerscfl found a more realistic value E g (T) = 7.11 cV and E g {X) = 9.11 cV through a Hartree- 
Fock calculation including correlation at the second order in perturbation theory. These results confirm the general 
findings in other semiconducting and insulating materials, that DFT-LDA underestimates the fundamental gap, while 
the reverse happens in Hartree-Fock calculations, although the inclusion of correlations in the latters improves the 
agreement with experiments. Our GW corrected bands at the experimental lattice parameter, yield gaps at X and T 
equal to E g (X) — 6.45 eV and E g (T) — 7.12 eV respectively, presenting a better agreement with experimental data. 



SrO 


a(A) 


B(Mbar) 


E coh (eV) 


Exp« 


riment 


5.16 


0.906 


10.45 


LDA 
LDA 


31 

57 


5.25 
5.1-5.18 
5.22 


1.06 
0.88-0.82 
1.07 


10.9-10.5 
9.59 


Present work 


5.07 


1.04 


11.9 



TABLE IV. Structural properties of Cubic SrO. A comparison is made with Hartree-Fock (Ref. |50[ ) and LDA calculations 
(Refs. |l]j37|). 
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SrO 


E e (T) 


E g (X) 


E g (L) 


Minimal Gap 


Experiment 
Ref. [i]] 
Ref. ~ 


5.896 
6.08 


6.28 
5.79 




5.9 (r-r) 


LDA (Refj37|) 
X a (Ref. |§) 
HF+corr. (Ref. Q 


4.26 
5.10 
7.11 


4.03 
9.11 


7.3 
12.36 


3.8(F-X) 
3.9(r-X) 
8.54(F - T) 


Present work 

LDA 1 " 

LDA* 

LDA+GW f 
LDA+GW* 


4.29 
3.92 
7.54 
7.12 


3.03 
3.11 
6.39 
6.45 


7.79 
7.47 
11.17 
10.81 


3.00 (r - x) 

3.04 (r - X) 
6.37 \v-X) 
6.39 (r - X) 



TABLE V. Gap of Cubic SrO (in eV), at the high-symmetry point of the Brillouin zone. The gaps are calculated either at 
the theoretical lattice parameter (') or at the experimental one (*). 
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VI. STRONTIUM TITANATE 



The computed ground-state properties of the cubic phase of SrTiC>3, such as the equilibrium lattice parameter^hii 
bulk modulus and the cohesive energy, are showed in Table VI. A good agreement with previous LDA calculationsEj'Ej 
is found. We stress that in all these calculations the selfconsistent electron density includes the semicore states (Sr 4s, 
Sr 4p, Ti 3s and Ti 3p). With respect to experimental data, a slight underestimate (-1.3%) of the lattice parameter 
is obtained, while the computed cohesive energy is ~ 20% larger than that measured, as it is usual in the LDA. As 
far as the GW calculation is .concerned, we adopt the experimental value eoo = 5.82 of the optical dielectric constant 
of the cubic phase of SrTiOsE-l Both LDA and GW bands (Fig.||) are computed at the LDA theoretical equilibrium 
lattice parameter a t h- 

Firstly, we discuss the positions of the deep lying states not shown in Fig|| A deep, almost non dispersive, energy 
band is found at ~ — 54eV below the top of the valence band (VBM) in the LDA calculations, which originates from 
Ti 3s states. The inclusion of QP corrections in our GW-,scheme pushes it downward, at -58.9 eV below the VBM. 
This is in very good agreement with XPS measurements^, which show a peak around -59 eV. 

Analogously, in the LDA calculation, the narrow bands originating mainly from Ti 3p states are located at -30.5 
eV below VBM, very close to a band with dominant Sr 4s character, at about -29.5 eV. GW corrections push them 
at p33.7 eV and -32.7 eV, respectively. The XPS spectra show a 2 eV wide asymmetric peak centered at ~ —34.5 
eVEl As found for the Ti 3s band, the inclusion of quasiparticle effects greatly improves the agreement with the 
experimental data. I— . 

The very broad structure in the XPS dataH between -20 eV and -14 eV originating form O 2s and Sr 4p states 
is fully consistent with the position of our GW bands, while the LDA bands are higher in energy by about 3 eV on 
average. 

The valence band portion shown in Fig.|| has a dominant O 2p character with non negligible contributions from Ti 
3d states, mainly in its lower part. It is less sensitive to GW corrections, which almost rigidly shift the bands towards 
lower energies. As a consequence, the LDA and GW bandwidths along the A line (4.15 eV and 4.11 eV, respectively) 
are both in very good agreement with the value of 4.2 eV measured in angle-resolved photoemissionH. Similarly, 
the total bandwidth Ri v — R±2 jV (see also Fig.||) turns out to be 4.8 eV in the LDA and 4.9 eV with quasiparticle 
corrections. 

A major difference between the LDA and GW results regards the direct band gaps, which roughta-differ by about 3 
eV. The onset for optical absorption was measured at 3.2 eV by CardonaE3 and at 3.34 eV by BlazeyES, and interpreted 
as a T\5 tV — > T25' , c transition. The strong absorption peaks at ~ 4 — 5 eV were attributed to vertical transitions at 
X and M, respectively. On the other hand, the fundamental gap derived from combined photocmission and inverse 
photoemission spectra, which is free from excitonic effects, was estimated to be (3.3 ±0.5) eVE3. With respect to these 
experimental values, the LDA E g (T) is underestimated by almost 1 eV, while the GW one is overestimated by about 
2 eV (see Table VI) In this respect, the inclusion of self-energy corrections at a perturbative level does not improve 
the LDA. This issue is discussed in the next section, in relation to the validity of our model screening for transition 
metal oxides such as SrTi03. 



SrTiOs 


a(A) 


B(Mbar) 


E coh (eV) 


Experiment 


3.903 


1.83 


31.7 


Ref. 


52] 


3.870 


1.94 




Ref. 


53 


3.864 


1.99 




This work 


3.850 


2.03 


37.88 



TABLE VI. Structural properties of cubic SrTiOs. Both Refs. ^2j,g3| use the LDA with ultra-soft pseudopotentials. The 
experimental values of a and E co h are quoted in Ref. pq and B is taken from Ref. B9J. 
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SrTiOa 


E e (T) 


E e (X) 


E g (M) 


E e (R) 


Minimal Gap 


LDA — 

LDAH 


2.16 
3.77 


4.1 




5.2 


1.79 (R-T) 
2.88 {R-T) 


This work 
LDA 

LDA+GW 


2.24 
5.42 


2.85 
6.10 


4.17 
7.29 


4.82 
8.15 


1.90 (R-r) 
5.07 (R-T) 



TABLE VII. Computed E g for direct transitions from the VB to the CB of cubic SrTiOs, in eV. The minimal gap is also 
given. 
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r x m r r x 

FIG. 3. Upper part of the valence bands and lower part of the conduction bands of cubic SrTiOs, at the theoretical 
equilibrium lattice parameter. Dashed lines: LDA bands. Solid line: GW bands. 



VII. DISCUSSION 



The agreement of our numerical results with the available experimental data is rather inhomogeneous: while for 
MgO and SrO a reasonable agreement is found, for SrTiOs the QP bands around the Fermi level seem to be worse 
than the LDA ones. It is the aim of this section to discuss our results with respect to two points: On the one hand, 
the quality of the ground-state calculation from which the perturbative GW calculation starts. On the other hand, 
the reliability of our model dielectric function used in the calculation of the QP spectra of the oxides considered here. 

We have found that the values of the direct gaps are very sensitive to the value of the lattice parameter. A good 
description of the ground state structural properties is thus a prerequisite to obtain QP spectra that can be compared 
to experimental data reliably. This is especially the case of MgO and SrO. In Table VIII the relative variations A g (fc) 
of the direct gaps, both in the LDA and in the GW approximation, are given as a function of the lattice parameter. 
One can note that: (i) The largest variations of the direct gaps A g (fc) are nearly always obtained at the LDA level. 
(ii) Ag(fc) depends on the actual k point, consistently with the character of the wavefunction. The largest A g (fc) are 
generally found for the strongest bonding-antibonding combinations across the gap. (Hi) Increasing A g (r) are found 
for SrTiOs, SrO and MgO, in ascending order. 

One should thus be careful when discussing the quality of the self-energy corrections to the electronic structure, 
since it may be strongly dependent on variables, such as the lattice parameter, which are external to the theory itself. 
Our data for MgO, obtained at the pap«rimental lattice constants a exp , compare well with other more sophisticated, 
ab-initio GW results also using a eX p"^, as far as the direct gaps are concerned (see Table HI). Both for MgO and 
SrO, a fair accordance with experimental data is found when adopting a exp . We thus conclude that, still at the level 
of the ground state structural properties, one should obtain a good agreement with experimental data (for example, 
by using generalized gradient corrections to the LDA) before starting the calculation of the quasiparticle energies. 

Care must also be taken when comparing the theoretical gaps and those derived from optical measurements, because 
of the existence of excitonic effects in the latters, not accounted for in the present theory. For MgO, these effects should 
be rather small, as suggested experimentally by Roessler and WalkeiCJ and calculated by Benedict and coworkers^. 
For SrO and SrTi03 there is still a lack of ab initio calculations of one- and two-particle excitations to address fully 
this issue. 

Another important point is the quality of the quasiparticle corrections, in relation to the model dielectric function 
e(q, p) used in the GW calculation. As it is evidenced by Eq.^, the screening depends on the value of the optical 
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dielectric constant of the oxide. As in the previous discussions regarding the lattice parameter, one can use 
either the experimental optical dielectric constant available in the literatuK®, or compute it in the framework of the 
self-consistent theory adopted, such as the DFT. In a recent paper, Shirleyc3 obtained for MgO a theoretical optical 
dielectric constant, at the experimental lattice parameter, e^ 1 ^ = 3.03 in the LDA, to be compared with the measured 
£oo = 2.95E3. Extrapolating our results - obtained with e^ = 2.95 and = 3.50 (see Sec. IV) - we can conclude 
that our quasiparticle energies would vary by less than 0.1 eV when passing from eoo = 2.95 to = 3.03. However, 
a more satisfactory theoretical framework would be to compute e m from first principles, in conjunction with the 
use of an approximation for the exchange-correlation energy, capable to reproduce the experimental lattice constant. 
Although in the present case the results would not be much affected, this would open the way to the calculation of 
quasiparticle excitations of systems for which no experimental measurements of the optical dielectric constant are 
available. On the other hand, the use of model dielectic functions might open the possibility to treat complex systems 
consisting of many inequivalent atoms, such as low-symmetry crystals, surfaces and heterostructures. 

Considering the numerical results obtained, and the previous discussion, we note that the model dielectric function 
works reasonably well for both MgO and SrO, which are compounds with fund amental gaps much larger than the 
semiconductors ZnSe, GaN, and SiC previously studied with the same methodcilO. We thus conclude that the validity 
of the screening model (EqJ^) is neither a function of the gap value nor of the ionocovalent character of the bonding 
in the crystal under study. On the other hand, it fails when applied to SrTiOs, a transition metal oxide characterized 
by a moderate optical gap (slightly larger than 3 eV according to Refs. |6j]56|). 

In order to clarify the reasons for such a failure, we computed the bare exchange contribution to the QP gap at the 
r point, for both MgO (10.46 eV) and SrTiO,3 (12.28 eV), in first-order pertubation theory. This difference stems 
from the different character of the states at the bottom of the conduction bands, which have a Ti 3c? character in the 
case of SrTiOs, and are much flatter than those of MgO, having a dominant Mg 3s character. This is fully consistent 
with the stronger localization of the Ti 3d states, which could make our treatment of local-field effects (see Eqs.|| 
and 0) rather inappropriate. In fact, the choice of the effective density p n ^ (Eq.||) takes into account the different 
localization of the states only through their superposition with the electron density p(r). Moreover, its definition 
contains an ambiguity, since its precise value depends on the core-valence partitioning of p{r). It is interesting to 
consider the extreme limit in which e(g, p) is no longer q dependent and is simply replaced by the optical dielectric 
constant £oo. In this limit of screening, the only contribution to the GW correction to the DFT-LDA gap comes 
from the screened-exchange term, since the Coulomb hole gives a rigid shift equal for all the bands. Interestingly, 
this correction provides gap values equal to 4.3 eV for MgO and 2.7 eV for SrTiOs, in much better agreement with 
experiments in the latter case than when using e(q, p n ^) given by EqJ|. This points out the need for a better inclusion 
of local field effects in the model dielectric function for transition metal oxides, such as SrTiOs. Another source of 
failure relatively to SrTiOs could be also the use of a perturbative GW scheme: a possible solution to this issue would 
be the use of an efficient self-consistent GW approach which has shown to give good results even in systems in which 
d orbitals play a fundamental roleS. 





MgO 


SrO 


SrTiOs 


(ath — a C xp)/a C xp(%) 


-2.0 


-1.7 


-1.3 


Ag UA ( r ) 


-6.4 


-5.4 


-1.7 


A? w (r) 


-4.1 


-3.6 




A g — (X) 


-0.7 


1.5 


-2.2 


A g GW (X) 


-0.7 


0.6 






-2.7 


-2.5 


-4.4 


A« W (L) 


-2.1 


-2.0 





TABLE VIII. Dependence of the direct gaps E s (k) of MgO, SrO and SrTiOs at high-symmetry points k of the Bril- 
louin zone as a function of the relative variations of the lattice parameter (ath — a C xp)/a C xp. The function A g is defined as 



8 ' , — g \ a " xp — — . The values are computed either by using the LDA, or by including quasiparticle corrections within 
the GW approximation. 
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VIII. CONCLUSIONS 



We have calculated the gound state properties of cubic oxides MgO, SrO and SrTiC>3 within the DFT-LDA, 
and applied a perturbative GW method, based on a dielectric model which includes approximately local field and 
dynamical effects, to determine their quasiparticle energies. We have tested the dependence of the calculated spectra 
on the parameters entering the calculation, such as the optical dielectric constant and the lattice parameter a. 
In particular, we have shown that quasiparticle energies in alkaline-earth oxides are sensitive functions of the lattice 
parameter, even at the LDA level. As a consequence, the use of either the experimental or the theoretical values for 
a can influence the quality of the results whenever a t h differs from a oxp . 

Through a careful comparison with the available experimental data and previous ab-initio GW calculations, we 
have showed that the simplified GW method works reasonably well for systems, such as MgO and SrO, in which 
the electronic states of the valence band and the bottom of the conduction band mainly consist of s- and p-like 
states. Moreover, in the three oxides, GW corrections have been found to be crucial to account for the energies 
of semicore electron states, such as O 2s, Ti 3s, Ti 3p, Sr 4s and Sr 4p, which agree well with the peak positions 
recorded in photoemission spectra. When localized d— like states contribute to bands around the Fermi level, our 
method produces an overestimation of the self-energy corrections. This is the case of SrTi03, for which it is not clear 
at this stage whether a higher order perturbative approach, in which the LDA wavefunctions are renormalized, is 
needed, or whether this failure calls for a more refined treatment of local-field effects. 

We thank Lucia Reining for fruitful discussions and careful reading of the manuscript, Tristan Albaret for his help in 
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for MgO. Computing resources were granted by the IDRIS/CNRS computational centre in Orsay, under project 
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